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We use lattice QCD simulations to test some of the predictions of proposed AdS/QCD (holo- 
graphic) duals for QCD. In particular, these duals predict that the scale of chiral symmetry break- 
ing (^SB) can be varied independently from that of confinement, with the proviso that the scale 
of ^SB cannot be longer than that of confinement. We simulate lattice QCD with 2 quarks in 
the fundamental representation of colour and with additional 4-fermion interactions (suggested 
by AdS/QCD), at finite temperatures. For sufficiently strong 4-fermion interactions, the decon- 
finement and ^SB transitions occur at different temperatures, the separation depending on the 
4-fermion coupling. This confirms that the scales of confinement and %SB are, in general, differ- 
ent. 
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1. Introduction 

With standard actions, simulations of the finite temperature phase structure of lattice QCD with 
fundamental quarks indicate that the deconfinement transition and the chiral-symmetry restoration 
transition are coincident \M. For quarks in other representations of SU (3) colour, deconfinement 
occurs at a lower temperature than chiral-symmetry restoration [Q, f3j. 

On the other, hand suggested holographic duals of QCD with fundamental quarks, inspired by 
AdS/CFT duality [0, 0, H], predict that the scales of confinement and of chiral symmetry breaking 
C^SB) can be varied independently [0, §]. However, since confinement produces ^SB, the length 
scale associated with #SB cannot be greater than that associated with confinement []9[ [RJ. 

These proposed holographic (string/gravity) duals of QCD suggest that the scales of confine- 
ment and ^SB can be decoupled by the addition of (non-local) 4-fermion interactions to QCD. 
Adding attractive 4-fermion interactions binds the quarks and antiquarks more tightly so that qq 
pairs can condense and spontaneously break chiral symmetry at shorter distances. 

We consider lattice QCD with 2-flavours of staggered quarks and local 4-fermion interactions 



of the Gross-Neveu Ql lQ/Nambu-Jona-Lasinio [ |12| , |13Q type. We work at finite temperature and use 
the deconfinement and chiral-symmetry restoration temperatures as our measure of the scales of 
confinement and #SB. These are measured as functions of the 4-fermion coupling. These simula- 



tions are described in reference [14] 



2. The lattice action 



The lattice quark action is the standard #QCD action []15P, which is the traditional staggered 
quark action augmented by a chiral 4-fermion term. Expressed in terms of the auxilliary fields a 
and % to render it quadratic in the fermion fields it is 



Nf/4 

: E E*/ 

f=l s 



£ = (— \y + y +l+t and 7 is the inverse 4-fermion coupling. This preserves the exact U{\) axial 
flavour symmetry of staggered fermions. 

Simulations are performed using the exact RHMC algorithm [[if]] to tune to 2 flavours. The 
deconfinement transition is determined as the position of the rapid increase in the Wilson Line 
(Polyakov Loop). The chiral-symmetry restoration phase transition occurs where the chiral-condensate 
YY vanishes (m = 0). 



3. Simulations and Results 

If we turn off the QCD interactions, we are left with a 4-fermion model which has a bulk 
transition at 7= 7. 1.7. At finite temperature, this will be shifted to smaller 7 (stronger coupling). 

Our finite temperature simulations are performed on N t = 4 lattices. We keep 7 > 7. so that at 
high temperatures - weak gauge coupling (large /3 = 6/g 2 ) - the theory is in the chiral-symmetry 
restored phase. Previous simulations at 7 = 10 and 7 = 20 indicate that the deconfinement and 
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chiral-symmetry restoring transitions are coincident [17]. Hence 7 > 10 represents a weak 4- 
fermion coupling. We simulate at 7 = 2.5, a strong 4-fermion coupling and 7 = 5, and intermediate 
coupling. Nf = 2 and m = 0. At 7 = 2.5 our lattice sizes are 16 3 x 4, 24 3 x 4 and 32 3 x 4. At 7 = 5 
our lattice sizes are 12 2 x 24 x 4, 24 3 x 4 and 32 3 x 4. Typical run lengths for each parameter set 
are 50,000 or 100,000 trajectories. 

Figure [l] shows the Wilson Line and chiral condensate as functions of j8 for the two 7s con- 
sidered. At 7 = 2.5 the deconfinement transition and the chiral-symmetry restoration transition are 
well separated. At 7 = 5 the two transitions are close, but still clearly separate. It is clear that the 
finite size effects are small near the deconfinement transition. As expected, the finite size effects 
are considerable close to the chiral-symmetry restoration transition. These graphs include points 
(and for 7 = 2.5 a lattice size) in addition to those in our earlier publication [p"4|]. 
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Figure 1: Wilson line and chiral condensate as functions of j3 for (a) J — 2.5 and (b) y = 5 in lattice units. 



At 7 = 2.5 we estimate the position of the deconfinement transition from the peak in the Wilson 
Line susceptibility, using Ferrenberg-Swendsen reweighting [18] from j8 = 5.545. This yields 
fid = 5.547(3). This is possible because the distributions of plaquette values from the simulated /3s 
in the neighbourhood of this transition overlap (see figure ||a). The j8s for the simulations near the 
chiral transition are not close enough for such an estimate (see figure |3|a). Our (more subjective) 
estimate for this chiral transition is (5 X =6.85(5). 

For 7 = 5, we estimate the position of the deconfinement transition to be jS^ = 5.420(4). For 
the chiral transition we estimate j3 x = 5.450(5). 

Figures HJ3] show the plaquette distributions near the deconfinement and chiral transitions. 

The deconfinement j8, /3j, is restricted to the range pd{7 = °°) < Pd < Pd(j = 7c) with j8^(7 = 
7c) < fid{quenched) . Hence 5.25 < fid < 5.6925. [The lower bound (zero 4-fermion coupling) 
is from Jl9[ [20] ] while the upper limit (quenched QCD) is from E2|].] The chiral-symmetry 
restoring j3, fi x , is in the range fid < fix — 00 > where the lower bound is achieved for small 4- 
fermion coupling, and the upper for y = y c . Table [j] summarises these results. (Note that we have 
not included error bars on the results of others.) 
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Figure 2: a) Plaquette distributions close to the deconfinement transition on a 24 3 x 4 lattice at 7= 2.5. 
(Since we have twice the statistics at j3 = 5.545 as at the other j3s, we have divided these points by 2). b) 
Same, but for 7=5. 
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Table 1: Deconfinement and chiral-symmetry restoration j3s as functions of 7. 



4. Discussion and Conclusions 

• Adding extra 4-fermion interactions with sufficient strength to the lattice QCD action sep- 
arates the deconfinement and chiral-symmetry restoration transitions at finite temperatures. 
We are able to change the separation of the scales of confinement and #SB by varying the 
4-fermion coupling, as predicted from proposed holographic duals of QCD. 

• 5.25 < p d < 5.6925, while fa < p x < 00. 

• At 7 = 2.5 - strong 4-fermion coupling - T x ~ 107^. At 7 = 5 - intermediate 4-fermion 
coupling - T x ~ 1.047^. For 7 > 10 - weak coupling - T x = Tj. 

• We should perform simulations at more 7 values in the range 2 < 7 < 7. Our simulations 
used N t = 4. A more complete study should include other Af, values. 
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32 J x4 lattice 7 = 2.5 
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Figure 3: Plaquette distributions close to the chiral-symmetry restoration transition on a 32 3 x 4 lattice (a) 
at y= 2.5, (b) at y= 5. 



• We have qualitative agreement with holographic QCD. We need to make the comparison 
more quantitative. 

• Is the deconfinement transition a phase transition, or merely a crossover? As the 4-fermion 
coupling increases to infinity, we expect the theory to approach quenched QCD, where de- 
confinement is a first order transition. We also see some indication of this transition be- 
coming stronger as the 4-fermion coupling weakens and the two transitions approach one 
another. 

• With a little more work we should be able to determine the universality class of the second 
order chiral transition. 
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• The two phase transitions appear to coalesce at a finite (non-zero) 4-fermion coupling (y a 
little larger than 5). 

• We used a local irrelevant 4-fermion interaction, which will not survive the continuum limit. 
Does the non-local 4-fermion interaction indicated by the AdS/QCD models survive the con- 
tinuum limit, i.e. does it define a non-perturbatively renormalizable theory? An intermediate 
model would be a Yukawa model where the auxilliary fields have full 4-dimensional scalar 
dynamics. Such a model is (or can be) perturbatively renormalizible, but is believed to be 
non-perturbatively trivial. 

• Theories with quarks in higher representations of the colour group, where the stronger QCD 
coupling separates the confinement and #SB scales even without the extra 4-fermion terms, 
should be studied. The Bielefeld group have studied adjoint quarks where the transitions are 
separated and are required by symmetry to both be phase transitions [f||. Colour sextet quarks 



are also of some interest, especially with regard to conformal technicolor models (23]. 
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